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ABSTRACT 

We  Investigate  a one-dlmenslonal  loop  model  for  the  evaporative 
cooling  of  the  coronal  flare  plasma.  The  Important  assumptions  are 
that  conductive  losses  dominate  radiative  cooling  and  that  the  evapora- 
tive velocities  are  small  compared  to  the  sound  speed.  We  calculate 
the  profile  and  evolution  of  the  temperature  and  verify  that  our  assump- 
tions are  accurate  for  plasma  parameters  typical  of  flare  regions. 

The  model  Is  In  agreement  with  soft  x-ray  observations  on  the  evolution 
of  flare  temperatures  and  emission  measures.  The  effect  of  evaporation 
Is  to  greatly  reduce  the  conductive  heat  flux  Into  the  chromosphere 
and  to  enhance  the  EUV  emission  from  the  coronal  flare  plasma. 


I . Introduction 

It  1 s becoming  widely  accepted  that  flares  evaporate  large  masses 
of  chromospheric  material  Into  the  corona.  The  most  compelling  arguments 
for  the  occurrence  of  this  process  come  from  soft  x-ray  observations. 
Measurements  of  the  emission  measure,  which  appear  to  be  well  over 

105°  cm  ^ for  some  flares  (Phillips  and  Neupert,  1973).  indicate  a 

T 16 

total  mass  of  hot  (T  2 10  K)  plasma  In  the  corona  as  large  as  10  g. 

Similar  mass  estimates  are  obtained  from  loop  prominence  systems  which 
are  known  to  occur  In  post-flare  regions  (Jefferies  and  Orrall,  1965). 

This  amount  of  plasma  Is  too  large  to  condense  from  the  corona  and  must 
somehow  be  extracted  from  the  chromosphere. 

Neupert  (1968)  was  perhaps  the  first  to  realize  that  corona  flare 
plasma  is  the  result  of  evaporation  from  the  chromosphere.  Similar 
ideas  were  expressed  subsequently  by  Hudson  and  Ohki  (1972)  and  by 
Sturrock  (1973)*  In  these  papers.  It  was  Implicitly  assumed  that 
evaporation  occurs  during  the  Initial  deposition  of  energy  Into  the 
chromosphere  by  means  of  high-energy  electron  streams.  However,  there 
appears  to  be  observational  evidence  for  evaporation  during  the  decay 
phase  of  flares.  The  emission  measure  is  often  seen  to  increase  while 
the  total  thermal  energy  does  not  increase  and  the  temperature  decreases 
(Kahler  et  al.,  1970;  Horan,  1971 )•  "Moving  boundary"  types  of  models 
have  been  proposed,  to  explain  these  results,  by  Strauss  and  Papagiannls 
(1971)  and  by  Zauman  and  Acton  (197*0.  These  models  Involve  the  assumption 
that  flare  heating  occurs  only  at  the  top  of  a loop  of  coronal  plasma, 
which  results  In  an  expanding  hot  region.  The  loop  Is  assumed  to  be 
essentially  Infinite  In  length  so  that  the  effect  of  the  heating  on 
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the  chromosphere  is  not  included.  However,  these  models  encounter 
certain  difficulties:  in  particular,  they  cannot  account  for  the  origin 
of  the  large  masses  required  by  the  soft  x-ray  observations.  Further- 
more, if  the  "moving  boundary"  models  were  correct,  flares  should  heat 
the  corona  well  before  they  heat  the  chromosphere,  whereas  observations 
indicate  that  the  Ho  emission  of  flares  begins  well  before  the  soft 
x-ray  maximum.  We  therefore  take  the  view  that  the  Increase  of  emission 
measure  during  the  decay  phase  of  a flare  indicates  that  evaporation 
occurs  not  only  during  the  impulsive  phase  of  a flare  but  also  during 
the  decay  phase. 

From  the  above  considerations,  it  appears  that  there  may  be  two 
evaporation  mechanisms  occurring  during  a solar  flare.  One  is  associated 
with  the  direct  heating  of  the  chromosphere  by  the  energy  released 
during  the  impulsive  phase  of  a flare,  the  energy  probably  being  trans- 
ferred by  an  electron  stream.  This  stage  is  likely  to  involve  large 
pressure  gradients  and  large  velocities,  quite  possibly  leading  to 
shock  waves.  It  seems  likely  that  a separate  evaporation  process  may 
occur  after  the  primary  energy  release  is  finished,  this  evaporation 
being  driven  by  the  large  conductive  heat  flux  from  the  hlgh-tempera- 
ture  flare  plasma  contained  in  magnetic  flux  tubes  above  the  chromosphere. 
In  this  stage,  evaporation  is  a cooling  process  competing  with  radiation. 
It  seems  likely  that  energy  loss  by  evaporation  dominates  early  in  the 
decay  phase,  and  energy  loss  by  radiation  dominates  later,  when  the 
density  of  flare  plasma  is  higher  and  the  temperature  la  lower.  More- 
over, it  seems  likely  that  this  evaporation  stage  would  Involve  only 
small  pressure  gradients  and  velocities.  It  is  this  later,  "gentle", 
evaporation  process  which  we  discuss  in  this  article. 
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II.  Model 


Since  the  magnetic  field  strength  in  typical  poet-flare  regions 
appears  to  be  large  (B  2 10  g,  Rust  and  Bar,  1973).  we  assume  that  it 
dominates  the  plasma  so  that  mass  flow  and  heat  transfer  are  both 
parallel  to  the  magnetic  field.  The  geometry  which  we  adopt  is  that 
of  a flare  loop  consisting  of  a narrow  current-free  flux  tube  with  the 
independent  variable  s measuring  distance  along  the  loop  (parallel  to  B) 
from  the  top  of  the  loop.  The  cross-sectional  area  of  the  flux  tube 
is  given  by  A(s)  and  will  depend  on  the  form  of  the  magnetic  field. 

The  field  is  assumed  to  be  that  of  a point -dipole  or  of  a line-dipole, 
but  we  also  calculate  the  case  of  a loop  of  constant  cross  section 
corresponding  approximately  to  the  familiar  plane-parallel  model. 

We  make  two  assumptions  on  the  state  of  the  plasma.  The  first  is 
that  radiation  losses  are  negligible.  We  later  check  our  assumption 
by  comparing  the  radiative  cooling  time  with  the  conductive  cooling 
time  to  be  derived  below. 

Our  second  assumption  is  that  the  velocities  generated  by  the 
evaporation  process  are  small  compared  with  the  speed  of  sound.  Since 
the  temperature  is  high  enough  that  gravitational  effects  are  negligible 
(l.e.,  the  scale  height  is  large  compared  with  the  size  of  the  flux 
tube),  this  assumption  is  equivalent  to  setting 


.1  ■.  1. La  point  it  may  seem  peredoxicel  that  we  can  calculate  the 
plasma  velocities  by  neglecting  the  pressure  gradients,  since  these 
gradients  create  the  velocities  In  the  first  place.  However,  one 

4 


finds  that,  if  the  gravitational  term  is  negligible,  the  fractional 
variation  of  pressure  along  the  flux  tube  is  of  order  M^,  where  the 
Mach  number  M is  given  by 

M ■ v/C  , (2.2) 

where  C is  the  speed  of  sound.  Hence  our  assumption  (2.1)  is  equivalent 
to  the  assumption  that  M « 1,  which  must  be  verified  a posteriori. 

Equation  (2.1)  shows  that  the  pressure  is  independent  of  s.  We 
assume  that  there  is  no  heat  input  into  the  loop,  and  we  ignore 
radiation.  Hence  the  total  energy  in  the  loop  is  constant,  which 
Implies  that  the  pressure  is  also  Independent  of  time: 


p(s,t)  ■ const. 


(2.3) 


On  using  equation  (2.2),  one  finds  that  the  energy  equation  may  be 

expressed  as 


_ i £ (M 

2 p \at 


+ V 


ds/ 


(2.4) 


On  combining  equation  (2.4)  with  the  continuity  equation 

U + ifi  - 0 • (2*5) 

we  obtain 


h [*&"  ' Kll)]  " 0 • (2,6) 

We  see  from  (2.6)  that 

A(fpv  - Kg)  - f(t)  . (2.7) 
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The  temperature  of  the  chromosphere  (the  base  of  our  model)  Is 
sufficiently  low  that  the  heat  flux  F,  given  by 


F 


(2.8) 


may  be  neglected.  Similarly,  the  density  there  is  so  high  that  the 
mass  velocity  v may  be  neglected.  These  boundary  conditions  thus  lead 
to  the  simple  relation 


v 


2 K 2T 
5 p 3s  * 


(2.9) 


which  may  alternatively  be  expressed  as 


(2.10) 


Equation  (2.10)  can  also  be  obtained  from  the  aymmeti^  of  the  loop  since 
if  the  loop  is  symmetric  about  the  top,  it  follows  that  F and  v must 
vanish  there. 

The  heat  equation  may  now  be  expressed  with  only  temperature  as 
dependent  variable.  On  replacing  v In  equation  (2.4)  by  the  expression 
(2.9) , and  Chen  eliminating  p In  favor  of  T by  means  of  the  equation 
of  state,  we  obtain 


52  21.  iL 

2 T at  A 3s 


(2.11) 


We  use  the  Spltcer  (1962)  formula  for  thermal  conductivity, 

K - or  T^2 


(2.12) 


where,  with  sufficient  accuracy  for  our  purposes,  we  may  take  a - 10~^. 

Equation  (2.11)  is  nonlinear,  but  It  Is  aawnable  to  solution  by 
separation  of  variables.  We  assume  that  T Is  of  the  form 
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(2.13) 


T(s,t)  - TM0(t)<p(s)  . 

For  convenience,  we  choose 

0(0)  - 1 , <p(0)  - 1 , (2.1U) 


so  that  T Is  the  maximum  temperature  of  the  plasma  (the  Initial 
temperature  at  the  top  of  the  loop).  We  may  substitute  (2.13)  Into 
equation  (2.11)  and  obtain 


1 _JL 

2 7/2  dt 


• «*•«) 


"V 


where  ^(s)  Is  related  to  cp(a)  by 


♦ “ <P 


,7/2 


(2.16) 


In  equation  (2.15),  the  variables  are  aeparated  so  that  we  may 

p 

set  each  side  equal  to  -k  (independent  of  both  t and  a).  The 


equation  for  0(t), 


7/2 


(2.17) 


may  be  Integrated  to  give 


e(t)  - (i  ♦ t/T  r2/7  , 


(2.18) 


where 


_5 £_ 


2QfTM7/2k2 


The  equation  for  f becomes 


(2.19) 


0.2  1 /dlf  + 1 dA  11  . 

d#2  7 ♦ 'de  / A di  di 


- 


(2.20) 
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We  see  from  equation  (2.14)  and  the  symmetry  of  the  flare  loop  that 
we  require  that 


♦ ■ 1,  ^ ■ 0 at  s - 0 . (2.21) 

However,  if  sb  is  the  value  of  s at  the  "base"  of  the  model  (l.e.,  at 
the  chromosphere),  the  temperature  must  drop  to  the  low  chromospheric 
value  at  s^.  For  our  purposes,  it  is  convenient  to  set 

*(»b)  - 0 . (2.22) 

The  position  of  the  chromosphere  will  actually  be  not  at  the  point 
where  t|r  - 0 (i.e.  T - 0),  but  where  T ms  10^  K,  but  this  difference 
is  small  compared  with  the  size  of  typical  flare  loops.  The  conditions 
(2.21)  and  (2.22)  can  be  met  only  for  a unique  value  of  k in  equation 
(2.20).  This  value  will  of  course  depend  implicitly  also  on  the  area 
function  A(s).  For  either  type  of  magnetic-field  tube  (that  produced 
by  a point  dipole  or  that  produced  by  a line  dipole),  we  shall  wish 
to  consider  a range  of  values  of  sb>  Rather  than  prescribe  &b  and 
then  determine  k as  an  eigenvalue  of  the  differential  equation  (2.20 ), 
it  is  clearly  simpler  to  consider  a range  of  values  of  k and  determine 
sb  for  each  value.  This  may  be  done  by  finding  a solution  of  the 
equation  with  the  initial  conditions  (2.21),  and  then  finding  the  point 
sb  at  which  f ■ 0. 
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III.  Results 


(a)  Temperature  Profile 

For  the  case  of  a loop  of  constant  cross-section,  equation  (2.20) 
can  be  Integrated  directly  yielding: 


? /l  - - li1  F(P(,).  sin  75°) 

>/3  + 1 - ^ 3 ' 


+ 2.3lA  E (0 («p) , sin  75°)  = |v/f  ks  » 


where 


0(<p)  * arc  cos 


v/3  - 1 + va 

,>/3  + 1 - ^ 


and  F and  E are  the  elliptic  Integrals  of  the  first  and  second  kind 
respectively.  The  constant  k can  now  be  expressed  In  terms  of  the 
position  s^  of  the  base.  Using  equations  (3*1)  and  (3.2),  we  find  that 
for  the  planar  case, 


k = 1.60/s, 


For  loops  of  variable  cross-section,  equation  (2.20)  has  been 
Integrated  numerically.  The  Important  parameters  that  characterize  the 
geometry  are  the  loop  length  s^  and  the  ratio  of  the  cross-sectional 
area  A(o)  at  the  top  of  the  loop  to  the  area  A(sb)  at  the  base.  We 
term  this  ratio  the  "compression  factor"  T: 


r - A(0)/A(s.  ) - B(s.)/B(o)  . 


Rust  and  Bar  (1973)  find  that  f < 30  In  the  flare  of  August  7,  1972. 
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For  a given  form  of  A(s),  we  may  solve  equation  (2.20 ) for  cp(s). 
In  Figure  3-1  we  plot  the  temperature  scale  height  H implied  by  our 
model , 


H(T) 


I 

T ds 


(3.5) 


as  a function  of  temperature  along  the  loop.  This  is  done  for  the 
case  of  a constant-area  loop  and  for  a flux  tube  due  to  a line-dipole 
source  with  a compression  factor  T * 30.  The  results  for  a point-dipole 
source  are  almost  identical  to  those  of  a line-dipole.  The  fact  that 
the  scale  height  H becomes  infinite  at  the  top  of  the  loop  simply 
reflects  the  fact  that  the  spatial  derivative  of  the  temperature  is 
zero  at  the  top  of  the  loop.  This  peculiarity  is  without  physical 
significance;  it  simply  means  that  the  spatial  variation  of  temperature 
is  determined  by  the  second  derivative  of  the  temperature  rather  than 
the  first  derivative. 

The  true  size  scale  over  which  T varies  is,  of  course,  less  than 

the  loop  length.  For  T < 10  1 , the  scale  height  is  small  compared  to 

-2  5/2 

the  loop  length,  H(T)  < 10  s^,  and  decreases  as  T in  both  the 

constant-area  and  the  variable-area  cases.  The  loop  would  therefore 
appear  to  an  observer  to  be  almost  isothermal  with  a very  thin  transi- 
tion region,  as  is  the  case  for  the  quiet  corona. 

For  comparison  we  also  plot  in  Figure  3.1  the  temperature  scale 
height  given  by  the  static  model  discussed  by  Antiochos  and  Sturrock 

(1976).  From  this  figure  we  find  that  at  low  temperatures  (T  < 10  * ) 

2/5 

the  temperature  profile  is  given  by:  T “ s for  evaporative  cooling, 

2/7 

and  T « a for  static  cooling.  Hence  the  heat  flux  F is  linearly 
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proportional  to  T In  the  evaporative  models,  whereas  F ~ const.  In  the 
static  models.  This  result  argues  strongly  In  favor  of  the  evaporative 
model 8 . 

For  temperatures  and  size  scales  typical  of  flares,  the  heat  flux 
Into  the  chromosphere  implied  by  the  static  models  Is  very  large,  In 

Q — 1 

excess  of  10  ergs  cm  s . It  is  difficult  for  static  models  to 
account  for  the  dissipation  of  such  a large  flux;  however,  In  our 
present  model  this  flux  Is  reduced  by  the  evaporation  process.  The 
heat  flux  into  the  chromosphere  Is  approximately  three  orders  of 
magnitude  lower  for  an  evaporative  model  than  for  a static  model. 


(b)  Cooling  Time  Scale 

We  can  evaluate  the  cooling  time  Tc  In  equation  (2.18)  by  using 
equation  (3. 3).  For  typical  loop  parameters, 


10'  K,  n 


IQ10*5  cm  ^ and  s. 


109'5  cm 


(3.6) 


equations  (2.19)  and  (3*3)  yield 


o 

T *s  10  sec. 
c 

The  radiative  cooling  time  Is  given  by 


T 


R 


1 -B..  . 

2 n2MT) 


(3.7) 


(3.8) 


where  A(T)  has  been  computed  by  Raymond  et  al.  (1976).  For  the  above 

•a 

parameters  (3.6),  we  find  that  TR  ms  10  sec;  hence,  our  assumption  that 
conduction  dominates  the  cooling  is  justified.  However,  as  the  tempera- 
ture decreases  the  ratio  tc/tr  increases  until  eventually  radiation 
dominates  and  evaporation  becomes  unimportant.  For  the  loop  described 
above,  we  expect  this  to  occur  at  T < 5 X 10^  K. 
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We  have  also  calculated  the  cooling  times  for  loops  due  to  either 
a point  dipole  or  line  dipole  source.  The  effect  of  the  geometry  is 
shown  In  Figure  3.2,  where  we  plot  the  ratio  of  the  cooling  times  for 
a loop  and  for  a planar  (constant  cross-section)  model  as  a function 
of  T.  As  in  the  static  case  (Antiochos  and  Sturrock,  1976),  the  effect 
of  the  compression  factor  is  to  increase  the  cooling  time;  however, 
we  still  expect  that  conduction  dominates  over  radiation  during  the 
initial  stages  if  T > 10^  K. 

The  other  main  assumption  in  our  model  is  that  the  evaporation 
velocities  are  small  compared  to  the  speed  of  sound.  For  the  planar 
case,  the  Mach  number  can  be  obtained  directly  from  equations  (2.9) 
and  (3*1): 

* . f . to-10-6  «1/2  (i  - »3/2)1/2  (i  ♦ £-r6/7  0.9) 

2/o 

The  maximum  value  of  M occurs  at  t ■ 0 and  <p  ■ (2/5)  .5*  for 

the  parameters  in  (3.6),  we  find  that  M < lo”1,  so  that  our  assumption 
is  valid.  From  the  form  of  the  time-dependence  of  M and  T,  we  find 
that  M « T^  for  fixed  s,  hence  our  approximation  improves  as  the 
loop  cools. 

(c)  Radiation 

Various  authors  have  attempted  to  calculate  the  emission  measure 
and  the  temperature  of  flare  regions  from  soft  x-ray  observations 
(e.g.  Kahler  at  al.,  1970;  Horan,  1971).  Craig  and  Brown  (1976) 
point  out  that  these  calculations  are  inherently  inaccurate,  and 
suggest  that  one  should  Instead  assume  a definite  model  for  the 
structure  of  a flare  region  and  then  test  this  model  by  comparing  the 
expected  radiation  with  observations.  We,  therefore,  calculate  the 
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evolution  of  the  soft  x-ray  emission  predicted  by  our  evaporative  loop 
model. 


p 


A 


7 

Tj|  = 2 X 10  K,  these  values  for  and  correspond  to  a soft  x-ray 
range  from  1-2  keV.  We  also  plot  the  time  dependence  of  the  temperature 
and  of  the  emission  measure  in  Figure  3«3*  Note  that  even  though  the 
temperature  is  continuously  decreasing,  l increases  initially  due  to  the 
increasing  emission  measure  and  has  a maximum  at  t/fc  “ 50.  For  larger 
values  of  r^  and  r^,  e.g.  r^  = 2 (E^  = 4 keV)  and  r?  ■ 4 (E?  - 8 keV), 
the  exponential  term  in  the  emisslvity  (3.11)  dominates,  and  L decreases 
monotonlcally.  We  find  that  these  results  are  insensitive  to  the 
compression  factor  T because  the  x-ray  emission  comes  mainly  from  the 
high- temperature  material  high  in  the  loop. 

The  qualitative  features  of  Figure  3*3  are  in  agreement  with  soft 
x-ray  observations.  However,  we  would  not  expect  to  find  detailed 
agreement,  since  we  have  calculated  the  emission  from  only  a single 
loop,  whereas  a flare  is  likely  to  consist  of  many  loops  at  different 
stages  in  their  evolution.  In  addition,  we  have  neglected  any  heating 
which  may  be  present  during  the  decay  phase  of  flares.  However,  our 
results  do  predict  a rise  in  soft  x-ray  flux  concurrent  with  a decrease 
in  temperature,  whereas  in  the  static  models  the  emission  must  decrease 
with  the  temperature. 

We  believe  that  a strong  observational  test  of  our  model  may  lie 
in  the  spectrum  of  the  radiation.  From  the  form  of  the  scale  height. 
Figure  3.1,  we  note  that  an  evaporative  loop  contains  relatively  more 
plasma  at  low  temperatures  (<  10^  K)  than  does  a static  one,  which 
implies  that  the  spectra  of  the  emitted  radiation  (especially  at  EUV 
wavelengths)  should  be  significantly  different  In  the  two  cases.  This 
Is  clearly  Illustrated  In  Figure  3*4  where  we  plot,  for  both  swdels. 
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the  total  radiative  loaaea  «(T)  due  to  plasma  at  different  temperatures 
along  the  loop,  l.e. 

«(T)  - n2A(T)  H(T)  , (3.15) 

where  we  use  Cox  and  Tucker's  (1969)  values  for  A(T),  and  H and  n are 
In  units  of  loop  length  sb,  and  density  at  the  top  of  the  loop  n(s  ■ 0), 

respectively.  The  temperature  at  the  top  of  the  loop  is  assumed  to  be, 

7 

Tm  - 10  K.  We  note  that  in  the  static  case  the  highest  temperature 
plasma  is  responsible  for  most  of  the  energy  loss  whereas,  in  the 
evaporative  loop,  plasma  at  ~ 5 X 10^  K radiates  most  strongly.  In  a 
later  paper  we  intend  to  calculate  the  spectra  fox  both  models  in  detail, 
in  particular  the  emission  in  various  EUV  and  x-ray  lines,  and  to  compare 
these  calculations  with  recent  Skylab  observations. 
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Figure  3.i 


Figure  Captions 

: The  temperature  scale  height  as  a function  of  temperature 

along  the  loop  for  three  models.  The  top  solid  line  refers 
to  an  evaporative  model  of  constant  cross-section.  The 
broken  line  refers  to  an  evaporative  model  for  a loop  due 
to  a line-dipole  source  with  T - 30.  The  bottom  solid 
line  refers  to  a static  model  of  constant  cross-section. 

Ratio  of  cooling  times  for  a variable  cross-section  model 
as  a function  of  T.  The  solid  line  refers  to  a loop  due  to 
a point-dipole  source  and  the  broken  line  to  a loop  due  to 
a line-dipole  source. 

The  evolution  of  the  soft  x-ray  emission  X between  1-2 
keV,  the  temperature  T,  and  the  emission  measure  (E.M. ) 
for  a loop  of  Initial  temperature,  T^  « 2 X 10^  K. 

The  radiative  losses  per  unit  logarithmic  temperature  Interval 
of  plasma  in  a loop  as  a function  of  temperature  along  the 
loop.  The  -top  line  refers  to  an  evaporative  model  and  the 
bottom  line  refers  to  a static  model. 
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